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The  estimator  is  evaluated  while  operating  in  white  Gaussian  noise 
which  is  in  general  correlated  between  channels,  and  curves  are  given  for 
the  density  function,  mean,  and  variance  of  the  estimator  for  various 
noise  assumptions.  The  estimator  is  shown  to  take  advantage  of  the 
processing  gain  of  large  time-bandwidth  product  signals  to  reduce  the 
variance  of  the  time  delay  estimate. 
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CHAPTER  1 


INTRODUCTION 

The  problem  of  time  delay  estimation  has  received  much  attention  in 
the  literature*-6.  In  its  simplest  terms,  the  problem  is  to  estimate 
the  time  difference  of  arrival  of  similar  signals  in  two  different 
channels.  In  general,  the  literature  can  be  divided  into  two  distinct 
categories  passive  and  echo  location. 

In  the  passive  mode,  the  receiver  "listens"  to  a  source  in  each  of 
two  receiver  channels,  and  estimates  the  time  difference  between  the  two 
channels.  In  this  case,  very  little  may  be  known  about  the  form  of  the 
signal.  In  echo  location,  a  signal  is  transmitted  in  the  channels,  and 
the  receiver  "listens”  for  reflections  in  the  channel.  Here  the  signal 
form  is  known  to  be  a  (possibly  distorted)  time  delayed  replica  of  the 
transmitted  signal,  and  the  task  is  to  measure  the  difference  of  the 
time  delays  in  the  channels. 

In  each  of  these  modes  of  operation  there  are  two  common  methods 
used  to  estimate  the  time  delay:  cross-correlation  methods,  and  phase 
comparison  methods.  The  first  of  these  methods  performs  a  cross- 
correlation  of  the  two  received  signals,  selecting  as  the  estimate  of 
the  time  delay  that  value  which  maximizes  the  magnitude  of  the  cross- 
correlation. 

In  the  phase  comparison  method,  the  analytic  signal  from  one 
channel  is  conjugated  and  multiplied  by  the  signal  from  the  other 
channel  and  the  phase  of  the  resultant  product  is  averaged  over  the  time 


duration  of  the  signal.  In  order  to  assign  a  meaningful  estimate  of  the 
time  delay  from  the  phase  information,  the  signal  must  be  assumed 
narrowband. 

In  the  literature  cited  above,  these  estimators  and  variants 
thereof  have  been  extensively  analyzed  under  various  hypotheses. 

This  thesis  presents  a  method  of  time  delay  estimation  which  was 
originally  developed  by  Ricker^.  This  method  is  inherently  an  echo 
location  estimator  of  time  delay  which  uses  a  comparison  of  the  phase  of 
the  matched  filter  for  each  channel  to  estimate  the  delay.  It  will  be 
seen  that,  under  the  proper  assume tions,  this  estimator  is  capable  of 
using  large  time-bandwidth  product  signals  and  still  giving  a  meaningful 
estimate  of  the  time  delay  from  the  phase  information. 

Other  advantages  of  this  method  include  the  fact  that  it  handles 
Doppler-shifted  channels  with  ease,  and  that  with  proper  signal  design, 
it  can  isolate  multiple  scatterers  in  the  channels  and  estimate  the  time 
difference  for  each  scatterer.  With  this  estimator,  one  may  utilize 
signals  with  large  time-bandwidth  products  that  resolve  well  both  in 
time  and  in  frequency  so  that  one  can  simultaneously  estimate  the  total 
propagation  delay,  the  time  stretch  (Doppler),  and  the  time  difference 
of  arrival. 

The  next  chapter  will  present  the  receiver  structure  and  the 
estimation  procedure.  Chapters  3,  4  and  5  will  evaluate  the  perfor¬ 
mance  of  the  receiver  operating  in  additive  white  Gaussian  noise,  and 
Chapter  6  will  summarize  and  suggest  further  work  in  this  area. 


CHAPTER  2 


RECEIVER  STRUCTURE 


This  chapter  will  discuss  the  method  in  which  the  received  signals 
are  processed,  and  will  present  notation  and  assumptions  that  will  be 
used  in  the  chapters  that  follow.  The  receiver  structure  is  shown  in 
Figure  2-1.  It  is  assumed  that  a  signal  f(t)e^u>ot,  was  transmitted 
and  that  the  received  signals,  r^(tj,)  and  r2(t^),  are  samples  of  two 
time-delayed,  time-stretched  replicas  of  the  transmitted  signal  plus 
additive  noise.  The  function  f(t)  is  known  as  the  complex  envelope  of 
the  signal,  and  uQ  is  equal  to  2irfQ,  where  f  is  the  carrier  frequency. 
The  received  signals,  ^  and  r 2  can  be  written  as  follows 

rl(ti)  *  ^  f[s(t1-T1)]ej“°S(ti  V  +  nJ(ti) 

r2(ti)  "  ^  f [s(t1-T2)]ejaio6(ti  V  +  n2(ti) 

It  is  assumed  both  signals  are  of  the  same  energy,  and  that  the  complex 
envelope  is  normalized  such  that 


N 


1 


k  -  1,2  (2-1) 

so  that  in  the  absence  of  noise,  the  energy  of  each  sampled  signal  is  E 


The  time-stretch  factor,  s,  is  assumed  to  be  the  same  in  each  channel, 
and  is  related  to  the  Doppler  shift,  4>jj,  as 


(2-2) 
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The  Cime  delays  for  each  channel  are  denoted  by  and  T2>  and  the  dif¬ 
ference  in  time  delays,  At  -  is  the  quantity  to  be  estimated. 

The  noises  n^(t^)  and  n2^Ci^  are  a88ume<*  to  he  white,  complex  proc¬ 
esses  with  circular  symmetry8>9,  and  with  variance  of  real  and  imaginary 
2 

parts  equal  to  o^.  The  real  and  imaginary  components  for  each  channel 
are  denoted  as 

nk(ti)  -  xk(ti)  +  jyk(ti)  k  -  1,2. 

The  received  signals  are  processed  according  to  Figure  2-1.  The 
processing  signal,  r(t^),  is  a  time-delayed,  time-stretched  replica  of 
the  transmitted  signal  with  hypothesis  time-delay  and  time-stretch 
values  of  t  and  s  obtained  through  some  previous  estimation  procedure. 
The  processing  signal  is  denoted  as  follows: 


The  processing  signal  is  assumed  normalized  to  unit  energy,  i.e, 


N  , 

l  r(t  ) 
i-1  1  l 


,  N 
2-  I 

i-1 


flsU^T)] 


is 


In  the  absence  of  noise,  the  output  of  erch  summation  in  Figure  2-1 


xk  “  l  ?(ti)rk(ti) 
i-1 


/F  l  f*[s(t  1-0]fI.(t.-Tlt)JeJ“oI8<ti-|k>-®<ti-;)1  (2-3) 

i-1  K 


»  e  ^ 
OaIaL 


A  A 

which,  when  viewed  as  a  function  of  T  and  s,  can  be  considered  as  a 
two-dimensional  correlation  between  the  received  and  processing  signals. 

The  estimates,  t  and  s,  are  often  obtained  by  processing  the 
received  signal  with  a  matched  filter.  The  matched  filter  concept  is 
based  on  the  examination  of  the  magnitude-squared  of  a  function 
X(t,s,t,s)  over  an  appropriate  grid  of  ?  and  s  values,  choosing  as  the 
estimate  of  the  time-delay  and  time-stretch  those  t,s  pairs  which 
maximize  the  magnitude  squared,  | x( T »s » t,s) |  .  Note  that  for  signals 
with  sufficient  resolution,  the  matched  filter  approach  is  capable  of 
resolving  multiple  reflections,  identifying  the  t  and  s  values  of  each 
reflector  as  distinct  maxima  on  the  t,s  grid. 

An  important  function  that  arises  out  of  the  discussion  of  matched- 
filters  is  the  ambiguity  function.  The  ambiguity  function  depends  only 
on  the  particular  signal  employed,  and  is  equivalent  to  the  matched  fil¬ 
ter  for  a  noiseless  scattering  channel  consisting  of  a  point  scatterer. 
The  accuracy  with  which  one  can  estimate  time-delay  and  -stretch 
depends  on  the  behavior  of  the  signal  ambiguity  function  near  the 
origin.  For  the  reader  unfamiliar  with  matched  filter  concepts,  rather 
thorough  treatments  are  given  in  Van  Trees  10,  and  in  Cook  and  Bernfeld1^ 
for  the  narrowband  case,  and  in  papers  by  Altes^,13  for  the  wideband 
formulation. 

The  time  delay  estimator  of  Figure  2-1  works  as  follows.  First, 
the  received  signal,  r(t),  representing  a  composite  of  the  signals  rj(t) 
and  r£(t)  (or  a  related  signal)  is  matched-filtered  over  an  appropri¬ 
ate  t-s  grid,  yielding  |x(t,s,t,s) |^.  This  function  is  examined  over 


7 


the  grid ,  with  those  peaks  that  are  above  a  pre-selected  threshold 
representing  substantial  reflections  in  the  channel  with  different  time- 
delay  and  -stretch  values  estimated  by  the  corresponding  t  and  s  values. 

Once  a  peak  has  been  located,  the  i  and  6  values  corresponding  to 
the  peak  are  used  to  process  the  received  signal  for  channels  1  and  2, 
yielding  x1(t1,s,t,s)  and  *  For  each  Peak»  the  co®Plex  numbers 

Xj(t1,s,t,s)  and  x*^*8***8)  are  multiplied  together,  yielding  (in  the 
absence  of  noise) 

Q  -  XjXj  -  Eo.(|£*[J(t1-i)]£[s(t1-i1>]eJ"o<*'s,tl)- 
i 

*  E0 • (Ilf*  f »C  t±-T) ] s( t±-T1) ] 

f(s(tj-;?)]f*[s(tj-T2)]e:)“o(8‘s)(ti’tj))ejh)o8(T2"Tl)  (2-4) 

If  the  time-stretch  estimate  is  assumed  correct  (i.e  s-s) ,  the 
first  line  of  (2-4)  reduces  to 

Q  -  Eo*(p*(s(ti-T)]f  [stt^tp]  )• 
CIf[s(ti-^))f*[s(t1-T2)])e^o8(T2-T1> 

-  Eo-(A1A*)ejto  (2-5) 


where 


^  ■  If  [s(t±-T>]f [s(t±-Tk>] , 


♦o  ’  V'VV- 


Hence,  it  is  seen  that  if  the  time-stretch  estimate  is  correct,  the  time 
difference  of  arrival  can  be  obtained  from  the  phase  of  0  as 


if  the  quantity  (Aj^)  is  real.  Note  that  if  “  tj»  then  A2  *  Aj  so 

* 

that  the  quantity  AjA2  is  real.  However,  with  arbitrary  At,  Aj  and  A2 
must  be  real  independent  of  one  another  for  most  signals  of  interest, 
and  it  is  desired  to  find  sufficent  conditions  for  which  this  occurs. 

To  do  this  it  is  more  convenient  to  utilize  the  continuous-time  domain. 


In  this  domain 


,  * 

A^  »  Jf  (st-st)f (st-ST^)dt, 


where  the  energy  conserving  factors  have  been  ignored,  as  they  are 
unimportant  in  the  discussion.  Defining 


g(u)  -  f(u-a)  <->  G(u>)  -  e"J“aF(o>) 


h(u)  *  f(u-b)  <->  H(u>)  «  e  J  DF(u>), 


Parseval's  relation. 


/  g  (u)h(u)du  ■  /  G  (ci>)H(u))dw 


can  be  used  to  obtain 


A k  -  /  G*(u))H(u)dtu  -  /  |F(w)  | 2  e*1“8(Tk"T)du. 

—00  —00 

If  | F(— to)  |  *  |F(u>)|,  this  reduces  to 
00 

*  2  /  | F (to )  j  ^ '  cos  us(T^-T)du,  (2-6) 

o 

which  is  a  real  quantity.  Hence,  a  sufficient  condition  for  A^  to  be 

real  has  been  established.  If  the  complex  envelope  has  a  symmetric 

power  spectrum,  |  F(— to)  |  »  |F(to)|  ,  then  A^  is  real.  Note  that  this 

is  dependent  on  the  proper  time-stretch  estimate  but  is  independent  of 

the  delay  estimate.  Throughout  this  thesis,  it  is  assumed  that  this 

condition  is  satisfied  so  that  A}  and  A2  will  be  real  quantities. 

The  performance  of  this  estimator  depends  on  two  factors.  First, 

the  additive  noise  inherent  in  the  channels  will  degrade  performance  by 

adding  unwanted  terms  to  (2-5),  thereby  corrupting  the  phase,  <t>  . 

o 

A  A 

Second,  the  estimates,  x  and  s  will  not  exactly  match  the  true  values. 

A 

If  the  estimate  t  is  reasonably  close  to  the  true  values  x^  and  x^, 
performance  will  not  be  adversely  affected,  since  the  only  effect  of 
this  is  that  the  factors  A^  and  Aj  in  (2-5)  will  become  smaller  relative 
to  the  additive  noise  terms  as  can  be  seen  from  the  cosine  term  in 

.  A 

(2-6).  If  the  estimate  s  is  incorrect,  however,  the  estimator  becomes 
biased  due  to  the  additional  exponential  that  appears  in  (2-4).  The  x 
and  8  estimates  obtained  from  the  matched  filter  will  obviously  become 
less  accurate  as  the  noise  level  increases. 


In  Chapters  3  and  4,  the  first  of  these  problems  is  addressed,  i.e. 

that  of  the  unwanted  noise  terms.  In  these  chapters  it  is  assumed  that 

the  estimates  of  t  and  s  are  correct.  In  Chapter  5,  the  effect  of  a 

* 

T-s  mismatch  is  discussed  using  a  Cramer-Rao  lower  bound  anproach.  To 
simplify  the  analysis  and  notation,  it  is  assumed  throughout  the  follow¬ 
ing  that  the  channels  contain  only  one  scatterer. 


CHAPTER  3 


UNCORRELATED  NOISE 

The  first  case  to  be  considered  is  that  of  the  receiver  operating 

in  white  Gaussian  noise  uncorrelated  between  channels.  The  complex 

noise  process  is  assumed  stationary,  ergodic,  zero  mean,  circularly 

2 

symmetric,  with  variance  of  real  and  imaginary  parts  equal  to  on«  The 
real  and  imaginary  components  for  each  channel  are  independent,  and 
are  denoted  as  follows : 

n i < c i )  -  + 

n2(ti)  *  x2(t1)  +  jy2(tjL) 
or  more  succinctly  as 

nki  “  xki  +  Jyki  k  “  1»2 

where  the  first  subscript  identifies  the  channel,  and  the  second  is  the 
time  index. 

The  white  noise  assumption  requires  that 

E{xki  xkj}  *  E{yki  ykj}  “  °n  *ij  k  *  l'2 


where  is  the  Kronecker  delta 


1,1-j 

0,i*j 


Circular  symmetry  requires  that 

E{xki  yki)  -  0  for  all  i;  k  -  1,2 


and  the  additional  assumption  of  noise  uncorrelated  between  channels 
requires  that 

E^xlix2j  ^  ™  E^yli  y2j  ^  m  E^xli  y2j^  m  E^x2i  y 1 j ^  ® 


for  all  i ,J .  (3-1) 


The  magnitude-squared  of  the  noise  process  has  expectation 


E(|atl|2>  -  E{xJi  +  £)  -  2»2  k-1,2. 

The  received  signal  for  each  channel  at  time  t^  is 

rk(ti)  "  ^  fts(trTk)]ej“oS(ti"Tk)  +  Vti> 


or 


ki 


/r 

o 


ki 


eJ“os(ti-Tk) 


\i 


where  f^  =  ffsCt^-r^)]  is  the  complex  envelope  of  the  received  signal. 
If  the  processing  signal  r(t^)  is  denoted  as  follows: 


r(t.) 


fi  e 


-jU0S(ti-T) 


A  A  A  A 

where  f^  =  f  [s(t^-T)j  represents  the  complex  envelope  of  the  processing 
signal ,  then  the  output  of  each  summation  block  in  Figure  2-1  can  be 
written 


„  _  y  ;  _  ,  r  f  -j ws( t ,-t) 

<k  "  [  ri  rki  "  [  fi  6  °  1 


/T  fkieJ“os(ti-Tk)  +  nki 


/r  f.f. . 

o  i  ki 


ejw0[(s-s)ti+ST-sTk]  +  n^^g-jwosCti-T) 


/r 


(if  s-s) 


(3-2) 


where  the  last  line  assumes  the  proper  time-stretch  estimate,  s  *  s. 


Since  the  Gaussian  noise  components  are  all  assumed  independent  of 
one  another,  the  explicit  and  implicit  phase  factors  in  the  second  term 


of  (3-2)  can  be  absorbed  into  the  noise  process  without  altering  the 
statistics  of  the  situation.  This  yields 

*k  ■  ^  •  |  Vkl  +  IlfjKi  k  -  1.2.  0-3) 

The  second  term,  being  the  sum  of  independent,  zero-mean  Gaussians  is 
itself  a  zero-mean  Gaussian,  which  shall  be  denoted  by  n^*  The  real  and 
imaginary  parts  of  this  term  have  variance 


o 


2 


n 


A 

since  I  |f^|  ■  1,  so  that  (3-3)  can  be  written  as 

*k  ■  ^  \  ‘3“»S<,'Tk>  +  \  k  -  1.2. 

whet.  ^  -  E 

The  output,  Q,  of  the  estimator  of  Figure  2-1  is  then 


X1X2  " 


vr  A^V^'V  +  n. 


/T0  A2e'j“os(T'T2>  +  n* 


-  E  A. A_  eJ“os(T2"Tl)  +  /T  A-n.  +  /z~  A.n!  +  n.n!  (3-4) 

o  l  2  021  012  12 

where,  a6  above,  the  complex  exponentials  have  been  absorbed  into  the 
noise  terms  without  altering  the  statistics. 

It  has  been  previously  assumed  that  s-s.  If  it  is  now  also  assumed 
that  the  time  difference  of  arrival,  At,  is  small  and  that  the  time  delay 
estimate,  x  is  accurate  such  that  t  ■  x^  (k  -  1,2),  then  f^  ■  f^, 
so  that  A^  *  Z  |?  |^  ■  1, 


To  see  what  is  meant  by  the  words  “small”  and  "accurate,"  one  needs  to 


examine  (2-6)  where  it  is  seen  that  if  cos  a>maxs(Tk-T)  ■  1,  then  will 

be  at  its  maximum,  which  is  1.  Here  o>  denotes  the  maximum  frequency 

max 

of  the  complex  envelope  (i.e.,  a>  equals  one-half  the  signal  band- 

max 

width).  In  this  case  (3-4)  reduces  to 


Q  «  E  eJ“os(T2_Tl)  +  /E“  n.  +  /E"  n*  *  n.n*  (3-5) 

o  O  1  O  L  I  L 

The  presence  of  the  last  term  in  (3-5)  renders  further  analysis 
intractable,  so  it  is  desired  to  quantify  conditions  under  which  this 
term  becomes  negligible.  The  expectation  of  the  magnitude-squared  of 
of  the  process  can  be  found  as 

E{ |n1ri* | 2}  -  E{[Re(n1n*)]2}  +  E{[lm(n1n*)] 2} 


4  a 


where  E{‘)  is  the  expectation  operator, 
noise  ratio  is  defined  as 


Thus,  if  the  input  signal-to- 


h 


Input  Signal  Energy 
Expected  Input  Noise  Energy 


N*E{ (n 


ki1 


> 


then  the  ratio  of  the  magnitude-squared  of  the  first  term  of  (3-5)  to 
that  of  the  second  or  third  term  has  expectation 


E'VKI>  2° 


«  Nh 


(3-6) 


while  the  ratio  of  the  magnitude  of  the  first  term  to  that  of  the  last 
term  has  expectation 


Here,  as  in  the  previous  chapter,  N  is  the  number  of  input  samples.  By 
comparing  (3-6)  and  (3-7),  it  is  seen  that  the  last  term  of  (3-5)  is  of 
second  order.  For  example,  if  there  is  a  OdB  input  SNR  (h=l)  and  there 
are  1000  input  samples ,  then  the  second  and  third  terms  are  30dB  below 
the  first,  while  the  last  term  is  60dB  below  the  first.  For  the  rest  of 
the  development ,  it  is  assumed  that  the  quantity  Nh  is  large  enough  so 
that  the  last  term  of  (3-5)  can  be  ignored.  It  will  be  seen  that  for 
values  of  Nh  as  small  as  10,  this  approximation  yields  variances  con¬ 
sistent  with  those  obtained  through  computer  simulations. 

It  is  seen  that  for  large  Nh,  the  output  SNR,  given  approximately 
by  (3-6),  is  N  times  as  large  as  the  input  SNR.  The  factor  N  can  thus 
be  considered  as  a  processing  gain.  For  a  receiver  sampling  at  a 
frequency  equal  to  the  bandwidth  of  the  signal,  N  is  equal  to  the 
time-bandwidth  product  of  the  signal,  so  that  the  receiver  structure  of 
Figure  2-1  is  seen  to  have  a  processing  gain  equal  to  the  time-bandwidth 
product  of  the  signal. 

Defining  n3  such  that 

n3  -  /T  (rij+n*), 

(3-5)  reduces  to 

Q  -  Eoeju)os(T2“Tl)  +  n3,  (3-8) 

where  n3  is  a  complex  Gaussian  process  with  variance  of  real  and 


imaginary  parts  equal  to 


Q  can  be  written  in  terms  of  its  real  and  imaginary  parts  as 

q  m  e  e^o  + 

*  Eo(cos$o  +  jsin$Q)  +  u+jv 
«  a  +  u  +  j(b+v) 


where 


»o  -  V<VTi> 


a  =  E  cosd>  ,  b  =  E  sin<t> 

O  O  O  0 


u  “  Re{n3>,  var(u)  -  a ^ 


v  *  var(v)  * 


so  that  our  estimate  of  4> q  ,  denoted  by  ,  is 


In  the  absence  of  noise  it  is  noticed  that  4>+4>0  as  desired. 

In  order  to  evaluate  the  performance  of  this  estimator,  it  is 
desired  to  find  the  probability  density  function  for  4> ,  denoted  by  f ( 4> )  • 
The  density  function,  f ( 4> )  would  be  more  precisely  notated  as  f($|<fr  ), 
and  might  be  more  properly  regarded  as  a  likelihood  function.^  Also, 
it  is  noted  that  the  estimator  to  be  derived  is  a  maximum  likelihood 
estimator* 1 Qf  $0  for  the  uncorrelated  noise  case. 

The  numerator  and  denominator  of  the  arctan  argument  in  (3-10) 
are  independent  and  have  Gaussian  densities  with  mean  b  and  a, 
respectively.  The  joint  density  for  numerator  and  denominator  is  thus 
the  product  of  the  two  marginal  Gaussian  densities,  and  to  find  the 
density  for  4> ,  one  may  convert  to  polar  coordinates  as  follows,  and 
integrate  over  0  <  r<  ®.  Letting  u+a  =  rcos4>  and  v+b  *  rsin<|>,  the  joint 


density  becomes 


f (u,v) 


2ira; 


exp 


"I  /  ,  <-  \ 

-T—  (u  +V  )  -*• 


2a! 


f  (r ,  4> ) 


2ira: 


“1  2  2 

r  exp  — r —  [(rcos$-a;  +  (rsin$-b;  ]  , 


2a„ 


'UUiX 


where  the  factor  r  is  the  Jacobian  of  the  transformation.  The  desired 
density  is  obtained  by  integrating  over  r: 

00 

f <4>)  *  /  f (r ,4>)dr 
o 


2xo 


1  -1  r  2.^2 

~T~  exP  — J 


2oI 


[a^+bz]  *  /  r  exp  — —  [r2  -  2r(acosij»  +  bsin<(>)]dr. 

o  lot 


Recalling  that  a  ■'  E  cos$  and  b  ■  E  sin$  ,  and  using  a  well  known 

O  O  0  0 

trigonometric  identify,  the  above  exnression  becomes 

2 


fU)  - 


1 


2xo 


2^expT^ 


-E 

°  ■  /  r  exp  - — 5— [rZ  -  2E  rcos($-4>  )]dr, 
2o‘  o  2o*  00 


or,  making  a  change  of  integration  variable, 

2  2 
E  -  -E  - 

f ( d> )  - - —  /  r  exp  — [r  -  2  cos(4>-<t>  )r  +  l]dr, 


2xo^  o 


2  a; 


t  2  2  2 

which  is  a  function  only  of  and  ( — )  .  Recall  that  a,  *  2E  0  ,  and 

n  '  fT  4  n  «  ’ 


o  n 


h  ■  E  /2Na  ,  so  that 
o  n 


2E2a 
o  n 


Nh. 


Here  again,  h  is  the  input  signal  to  noise  ratio,  and  N  is  the  number  of 
input  samples.  Finally,  then,  the  density  is  written 


Nh=100 


Figure  3-2 


tm\ 


E{At} 


US  U  E 

o  o 


var{Ar}  ■  — - — 7r  c2 

(wj)2 

o 


-  .Ym  «V»  -  * 


where  the  first  integral  is  zero  due  to  its  odd  integrand  integrated 
over  symmetric  limits,  and  the  second  integral  is  1  since  it  represents 
the  area  under  density,  £(♦).  Hence,  with  this  choice  of  integration 
limits,  the  mean  is  equal  to  the  true  value  (i.e. ,  the  estimator  is 


unbiased),  so  that  the  variance  is  a  measure  of  the  fluctuation  of  $ 
about  the  true  value, 

The  estimate  of  time  delay  is  given  by 


which  has  mean  equal  to  the  true  value  (if  sas),  and  variance 


where  a  «  var(<t>)  is  shown  in  Figure  3-3. 


From  Figure  3-3,  it  is  seen  that  the  variance  of  this  estimator 


v,Y.  •• 

Y  v  V 

t'V.C 


decreases  approximately  as  1/Nh.  Under  previously  stated  assumption,  N 
can  be  replaced  by  the  time-bandwidth  product  of  the  signal,  and  the 
time  delay  estimator  is  seen  to  take  advantage  of  the  processing  gain  of 


k .  •  •  • .] 
-  •  .  «  1 


large  time-bandwidth  product  signals  to  reduce  the  variance  of  the  time 
difference  estimate. 


■m 

a  *  .  w  •  ' 

ISiF'l 


In  order  to  test  the  validity  of  the  variance  curve  of  Figure  3-3, 
the  estimation  procedure  of  Figure  2-1  is  simulated  on  a  computer.  The 
signal  used  for  the  simulation  is  based  on  a  Welsh  construction*?  with  a 
real  envelope.  The  signal  is  of  the  following  form 

f(  t)e^2irfot 


where 


f  - 

0 

30,000Hz 

i 

fi 

1 

550 

f(t) 

-  w(t)cos2irfit  (i-l)Ts  <  t  <  iTg 

2 

3 

750 

650 

TS  - 

subpulse  length 

4 

700 

.05  seconds 

5 

350 

6 

850 

w(t) 

-  50  dB  Taylor  Window*8 

7 

600 

8 

400 

9 

500 

10 

450 

11 

800 

12 

300 

The  ambiguity  function  of  this  signal  is  shown  in  Figure  3-4.  Note 
that  this  signal  has  good  resolution  both  in  T  and  in  s,  due  to  its 
large  time-bandwidth  product.  Note  also  that  by  giving  the  signal  a 
real  envelope,  one  obtains  a  symmetric  spectrum  as  required  by  (2-6). 

The  signal  is  sampled  at  1700Hz  and  matched-filtered  using  the 
correct  t  amd  s  values,  t»Tj»1.0,s"s*1.  Next,  1000  complex 


Gaussian  noise  samples  were  generated,  normalized  as  in  (3-9),  and  added 


to  0  as  in  (3-8)  to  obtain  a  sample  mean  and  variance  to  compare  with 
the  theoretical  unbiased  mean  and  the  variance  of  Figure  3-3.  The 


results  are  summarized  in  Table  3-1  table  for  $  ■  0  and  $  *  45°. 

o  o 

Table  3-1 

Theoretical  versus  Sample  Means  and  Variances 


I  Nh 

10 

100 

1000 

10.000  | 

Theoretical  Mean  (Deg) 

0.000 

0.000 

0.000 

0.000  I 

|  ♦  -  0* 
o 

Sample  Mean 

0.068 

0.045 

0.015 

0.005  1 

Theoretical  Variance 

380.0 

33.2 

3.3 

0.36  | 

Sample  Variance 

369.0 

34.0 

3.4 

0.34  | 

Theoretical  Mean 

45.00 

45.00 

45.00 

45.00  | 

1  *  -  45° 

1  o 

Sample  Mean 

45.84 

45.25 

45.09 

45.05  I 

l 

Theoretical  Variance 

380.0 

33.2 

3.3 

0.36  | 

1 

Sample  Variance 

379.6 

35.0 

3.5 

_ 

0.34  | 

It  is  seen  that  the  expressions  derived  above  for  the  density  func¬ 
tion  of  the  estimator  generally  predict  variances  consistent  with  the 
sample  variances  of  the  simulations  even  for  values  of  Nh  as  low  as  10, 
and  the  sample  mean  approaches  the  theoretical  mean,  for  large  values 
of  Nh. 

This  concludes  the  analysis  of  the  estimator  operating  in  uncorre¬ 
lated  noise.  The  next  chapter  will  discuss  the  considerably  more  com¬ 
plicated  case  of  noise  which  is  correlated  between  channels. 


CHAPTER  4 


CORRELATED  NOISE 


In  this  chapter,  the  receiver  performance  is  evaluated  when 
operating  in  white  Gaussian  noise  which  is  correlated  between  channels. 
As  in  the  previous  chapter,  the  noise  is  assumed  to  be  a  white,  cir¬ 
cularly  complex  Gaussian  process,  so  that 

E^xkiyki^  “  0  for  i5  k  “  1,2 

E{xkixkj}  "  E{ykiykj}  “  °n  6ij  k  "  1*2* 

The  noise  process  is  assumed  to  be  correlated  as  follows: 6 

E{xUx2J>  *  E(yU>V  *  pon  Slj 

E{xllV  ‘  -E{x2Jyll>  ■  lal  V 

A  physical  interpretation  of  p  and  A  is  given  in  Appendix  A.  It  is 
noted  here  for  future  reference  that  p  and  X  are  such  that 


P2  + 


A2  < 


1. 


(4-1) 


This  can  be  seen  from  examining  eigenvalues  of  the  covariance  matrix 

for  the  density  f (xj  ,yj  ^  ,  noting  that  since  the  covariance  matrix 

19 

is  positive  definite,  its  eigenvalues  must  all  be  positive. 

The  development  of  Chapter  3  may  be  followed  up  to  (3-2)  without 
modification.  Equation  (3-2)  is  repeated  here  for  convenience: 


X,. 


I  f.f„.  -  l  n,  . f ,  e-J“os<tl-,) 


1.2. 
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X  represents  the  output  of  the  k  b  summation  block  in  Figure  2— 1 ,  and 
(4-2)  assumes  the  correct  time-stretch  estimate.  Since  the  complex 
noise  is  assumed  correlated  between  channels,  more  care  must  be  taken 
when  absorbing  the  complex  exponential  of  the  last  term  of  (4-2)  into 
the  noise  process. 

If  the  second  term  of  (4-2)  is  written  in  terms  of  its  real  and 
imaginary  components, 


I 

i 


-jii)  S(t.-T) 

e  o  l 


uk + 


then  the  correlation  coefficients  for  u^ ,  ,  u2,  v2  can  be  expressed 

in  terms  of  the  input  correlation  coefficients,  P  and  X.  Using  the 
results  of  Appendix  B  and  the  fact  that  the  time  samples  of  the  noise 
are  independent,  it  can  be  shown  that  Uj ,  v} ,  u2,  v2  have  the 
same  correlation  coefficients  as  x^,  y^,  x2i,  and  y2J,;  that  is. 


E(u,u  }  -  E{v  v  }  »  po^ 

12  12  n 

2 

E{u.v0}  -  -E{u,v.}  *  Xo  . 
12  2  1  n 


It  is  noted  that  the  variances  also  remain  unchanged,  due  to  the  fact 
that  the  processing  signal  is  of  unit  energy. 

As  in  the  previous  chapter,  this  processed  noise  term  is  denoted 


by  V 


■  u^  +  ,1vk  *  I  nki  e  -^“o8^!  ^  k  •  1 ,2 


with  variance  of  real  and  imaginary  parts  equal  to  o^. 


If  is  again  defined  as 


l  fifki 


k  -  1,2, 


then  the  output  of  each  summation  block  in  Figure  2-1  can  be  written  as 


JT  a.  ,)V(K>  + 


xk  “  /Eo  \  ^  °'  k'  +  \  k  -  1,2 


so  that  the  output  0,  is  given  by 


Q  -  XjX*  ■  [■'iTAj  «JV(?-V  +  ^H/TAj  e-iV<;-V  +  ,*2] 


-  E  •A,A,ej"o:<T2-Tl)  A  /Ta,-,,.'3  A  /EVse3  V^'V 

ol2  o21  o  1  2 


A  n,n2 


(4-3) 


Defining 


*  n^e  jV(T  T2^  •  u^  +  jvj 


n2  “ 


_  — 1  to  s(t-t,  )  ”  .  - 

•  o  1  ■  U2  +  jv* 


it  is  again  desired  to  find  the  correlation  coefficients  of  u^ ,  v^, 
u„ ,  v„  in  terms  of  P  and  X.  The  results  follow  from  direct  application 


of  Appendix  B,  and  are  summarized  below: 


where 


E{u.u0}  -  E{v. v  }  -  (cp  -  dX)p2  -  po2 
12  1  2  n  n 

E{u  v  }  -  -E{v -U- }  ■  (cX  +  dp)o2  =»  Xo2 
1  Z  1  Z  n  n 


c 

d 


COS  U>  s(  T--T,  ) 

0  Z  1 

sin  u)  s(  )  . 
o  z  \ 


and  the  new  correlation  coefficients,  p  and  X  are  introduced  to  simplify 
notation  in  the  following  development.  It  is  noted  for  future  reference 
that 


P2  +  X2  -  (c2  +  d2)(p2  +  X2) 


2  2 
p  +  r  <  l 


where  p2  +  X2  <  1  as  in  (4-1). 


Equation  (4-3)  can  now  be  written  in  terms  of  and 


hj  as 


E  •A.A,e;,‘V,(T2"V  +  /e~ A~n, 

o  1  Z  o  2  1 


+  /eTa,  n*  + 
o 


12 


~  -*Us(T-h) 

hjhje  o  21 


■  Eo*e^u’o^T2  Tl^  +  /E^hj  +  Hj)  +  e^“os^T2~Tl^  ,  (4-4) 


where  the  last  step  assumes  that  the  estimate  of  time  delay  is  correct, 
and  that  the  time  difference  of  arrival  is  small  so  that  A^  *  1  and 
A.  ■  1,  as  in  the  previous  chapter. 


n.TrjxiV. *vi «n .vj^'yvi 


Again,  as  a  concession  to  tractability,  conditions  under  which  the 

last  term  is  negligible  are  quantified.  The  expectation  of  the 

•  -* 

magnitude-squared  of  the  process  n^  can  be  found  directly  through 
tedious  integration  to  be 


E{|^f*|2}  -  EdReC^*)]2}  +  E{tIm(^^)]Z} 

4(1  +  p2  +  X2)  <  8o4 
n 

-2  -2 

Where  the  inequality  arises  from  the  fact  that  p  +  X  <  1. 
with  the  input  signal-to-noise  ratio  defined  as 


z  :*m2i 


2*\ i 2i 


the  ratio  of  the  magnitude-squared  of  the  first  term  of  (4-4)  to  that  of 
the  second  or  third  term  has  expectation 


E((/r;k)2) 


(4-5) 


as  in  (3-6),  while  the  ratio  of  the  magnitude-squared  of  the  first  term 
to  the  expected  magnitude-squared  of  the  last  term  is 


EilnjSjl2)  8o‘ 


£  (-i-)2  -  £«*>2 

O 

n 


(4-6) 


From  (4-5)  and  (*-6)  it  is  seen  that  the  last  term  of  (4-4)  is  again  of 
second  order,  becoming  negligible  for  large  values  of  Nh.  For  example, 
if  Nh  -  1000,  then  the  second  and  third  terms  are  30dB  below  the  first, 


while  the  last  terra  is  at  least  57dB  below  the  first.  Throughout  the 
rest  of  the  chapter,  the  last  terra  of  (4-4)  is  dropped. 
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The  estimate  of  4»0  is  then  given  by  4> ,  where 


.  b  +  /E- (v,-v.) 

~1  r  o  1  2  _  -1  ^b  +  y ^ 

“  tan  [ - — - — ~  J  ”  tan  (— — 

a  +  /T(u,+u,) 


o'  1  2 

and  where  x  5  u^+  u2»  ^  £  vi“  v2*  In  order  to  evaluate  the 
statistics  of  this  estimator,  it  is  necessary  to  find  the  covariance 
matrix  for  the  numerator  and  denominator  of  the  argument  of  the  arc¬ 


tangent.  Denoting  this  by  R  ,  it  is  seen  that 
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It  is  noted  that  p  <1  since  p  +  X  <1.  The  joint  Gaussian  density 

*y 


is  then 


8 


Once  again,  the  change  of  variables  x  +  a  =  r  cos<(>,  y  +  b  *  r  sin<|>  is 
made  to  obtain  the  density  for  <j>  *  tan  After  considerable 

algebraic  manipulation  the  joint  density  f(r,4>)  is  obtained: 


where 
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The  desired  density,  f(<J>)  is  again  obtained  by  integrating  the  joint 
density,  f(r,<j>)  over  all  values  of  r,  i.e. 


f(4>)  -  /  f(r,<t>)dr. 
o 


This  integral  is  of  the  same  form  as  that  of  (3-11),  but  with  consid 
erably  more  complicated  coefficients  in  the  exponential.  These  coef 
ficients  can  be  substantially  simplified  by  writing  o  ,  c  and  p  1 
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Is  replaced  by  the  quantity  2Nh  to  get  the  density  in  terms  of  the  input 
signal-to-noise  ratio. 

2 

Roughly  speaking,  the  quantity  E0/an  controls  the  sharpness  of  the 
density,  while  the  correlation  coefficients  alter  its  shape.  The 
density  is  integrated  numerically  for  various  combinations  of  parameters 
yielding  Figures  4-1,  4,  7,  and  10. 

The  densities  are  then  used  to  generate  curves  of  the  mean  and 

variance  of  the  estimator  for  these  combinations  of  parameters,  yielding 

the  continuous  curves  in  Figures  4-2,  5,  8,  11  and  4-3,  6,  9,  12.  It  is 

again  noted  that  the  moment  quantities  depend  on  the  integration  limits 

chosen.  For  the  correlated  noise  case,  the  choice  of  these  limits  does 

not  seem  as  clear  as  it  was  for  the  uncorrelated  noise  case,  because  the 

estimator  is  now  biased.  If  the  same  limits  are  chosen  for  this  case 

however,  the  estimator  will  become  unbiased  for  large  values  of  Nh,  so 

the  integration  limits  are  again  chosen  as  $  -  ir  <  $  <  4  +  tt. 

o  o 

These  curves  are  again  verified  using  the  same  signal  and  sampling 
parameters  as  in  the  previous  chapter,  yielding  the  discrete  points  in 
the  mean  and  variance  curves.  The  noise  model  used  in  these  simulations 
is  described  in  Appendix  A,  with  the  incoherent  noise  correlation 
coefficient,  set  equal  to  zero  so  that  the  incoherent  noise  is  uncor¬ 

related  between  channels.  The  coherent  noise  is  the  same  in  each  channel 
except  for  a  phase  delay. 


-1d> 

e  c  n 


c 

li 


28  ■ 


of  E8tlmatot» 


OOI-MN 


Nh-100 


Figure  4-11 


Theoretical 


Variance  of  Estimator 


where  ■  tan  [— ) ,  p  and  X  being  the  input  noise  correlation  coeffi- 
cents.  The  incoherent  noise  has  variance  of  real  and  imaginary  parts 

a,2  -  <l-r)an2, 

while  the  coherent  noise  has  variance 


2  „  2 
a  ■  ra  , 
c  n 


/2  2  2 

where  r  ■  /p  +X  ,  and  is  the  total  noise  variance. 

From  these  figures,  it  is  seen  that  both  the  sample  and  theoretical 
means  approach  the  true  values  for  large  values  of  Nh.  Also,  the  vari¬ 
ances  predicted  by  theory  are  consistent  with  the  sample  variances  gen¬ 
erated  by  the  simulation  for  large  Nh.  The  theoretical  curves  are 
accurate  only  for  Nh  of  about  100  or  greater  due  to  the  approximation 
made  by  developing  the  last  term  of  4-4. 

Again  it  is  seen  that  for  large  Nh,  the  variance  drops  of  as  1/Nh. 
Replacing  N  with  the  time-bandwidth  product  of  the  signal  (as  in  the 
previous  chapter),  it  is  again  seen  that  the  estimator  takes  advantage 
of  the  processing  gain  of  large  time-bandwidth  signals  to  reduce  the 
estimator  variance. 

The  analysis  of  this  and  the  previous  chapters  has  assumed  that  the 


time-stretch  and  -delay  estimates  are  correct.  In  most  cases  these 
estimates  will  not  be  exactly  correct.  The  effect  of  a  T-s  mismatch  is 
the  subject  of  the  following  chapter. 


The  output  of  Figure  2-1  is  then  written  as 


n  >  v  V*  .  r  j ( 8—6 )  _  {  ,  jo  -  j*  *  -ja 

Q  xix2  o  e  0  e  r  fifii  6  1  :  fjf2j  e  j 
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(5-2) 


where 


>2  '  V<W  ' 


J* 


In  this  equation,  eJ  o  represents  the  true  signal  phase  difference,  and 


je 


it  is  seen  that  the  bias  represented  by  the  factor  eJ  cancels  itself. 


so  that  the  only  bias  remaining  in  the  first  term  is  that  contributed  by 


±10  ±1a 

the  double  summation.  It  is  also  noted  that  the  factors  e  and  e  i 


do  not  affect  the  correlation  coefficients  between  the  second  and  third 
terms  (see  Appendix  B).  The  last  term  will  be  dropped  under  previously 
stated  SNR  assumptions. 

If  the  assumption  is  again  made  that  the  time  difference  of 
arrival,  T2”T^  1®  ®wll  relative  to  the  bandwidth  of  the  signal  as  in 
previous  chapters,  then  f2i  -  f^,  so  that  the  first  term  in  (5-2) 


becomes 


VeJ*°  E  *ifii  «3°i  '  1  Kf2A  e~jaj  "  E*ej(W  |Z  f  f  ejai|2, 
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so  Chat  there  is  no  bias  in  the  estimate  of  $  due  to  an  incorrect  time- 

o 

stretch  hypothesis.  There  still  remains,  however,  a  bias  in  the  time 
difference  estimate.  This  bias  occurs  from  using  the  estimated  time 
stretch,  s,  rather  than  the  true  time-stretch,  s.  The  estimated  time 
difference  (in  the  absence  of  noise)  is 

AT  -  — is. - - t_  (1  -  £=1)  .  at(  1  -  szk)  (5-4) 

oi  a  o)  s  +  u  (s-;)  v  8 

o  o  o 

A 

where  At  is  the  estimate  of  t*ie  true  value,  4>q  »  * 

(i>os(t.-Tj),  and  where  a  binomial  expansion  yields  the  approximation. 

For  reasonably  good  time-stretch  estimates  (i.e. ,  |s-s|«l)  this  bias 
will  be  negligible. 

Comparing  the  magnitude-squared  factor  in  (5-3)  with  (2-3)  and 
noting  the  discussion  following  (2-3)  it  is  seen  that  this  factor  is  the 
magnitude-squared  of  the  matched  filter  output  whose  main  peak  is  at 

*  *  A  A 

t  *  Tj ,  8  ■  s,  evaluated  at  the  point  (t,s)  corresponding  to  the  esti- 
mates  of  the  the  true  values.  If  the  T  and  s  estimates  are  not  exactly 
correct,  the  magnitude-squared  factor  in  (5-3)  will  become  smaller, 
effectively  lowering  the  signal-to-noise  ratio. 

Thus,  It  is  seen  that  a  t-s  mismatch  does  not  affect  the  noise  terms 
(or  change  any  bias  due  to  the  noise  terms),  and  if  the  time  difference 
of  arrival  is  small,  the  only  effect  is  a  slight  bias  represented  by 
(5-4),  and  a  decrease  in  the  signal-to-noise  ratio  due  to  the  magnitude- 
squared  factor  in  (5-3). 


If  the  time  difference  of  arrival  cannot  be  assumed  small,  as  was 
necessary  to  derive  (5-3),  then  the  sum  must  remain  as  in  the  first  term 
of  (5-2).  The  problem  thus  becomes  more  difficult,  and  the  performance 


depends  on  the  particular  signal  employed. 

A  special  case  of  interest  here  is  that  of  a  real  envelope  signal 
(such  as  that  employed  in  the  examples  of  the  previous  two  chapters).  In 
this  case,  the  double  sum  in  first  term  of  (5-2)  can  be  written  as 
follows : 


l  fifli 


i  £  f*f* 
j  j  23 


e-J°j  -  £  £  f\f„  ej%(s-s)(ti-tj) 

i  j  J  ^  (5-5) 


(where  the  complex  conjugate  symbols  have  been  retained  for  future 
reference)  so  that  a  crude  upper  bound  can  be  placed  on  the  phase  bias 
contributed  by  this  term: 


l*bl  <  “0!s-s|T  (5-6) 

where  4>b  denotes  the  phase  bias  and  T  is  the  signal  duration. 

It  is  desired  to  obtain  an  idea  of  the  typical  size  of  the  estima¬ 
tion  errors  for  the  time-stretch.  The  estimation  errors  for  matched 
filter  receivers  have  been  considered  by  VanTrees10  and  by  Cook  and 
Bernfeld11,  who  use  a  Cramer-Rao  lower  bound  approach  to  derive  an 
expression  for  the  minimum  error  variances  for  this  estimate,  and  who 
show  that  the  estimation  errors  approach  the  Cramer-Rao  bound  for 
maximum-likelihood  estimators. 

In  general,  the  estimation  error  for  s  will  be  coupled  with  that  for 
t,  but  a  sufficient  condition  under  which  these  estimates  become 
uncoupled  is  that  the  signal  employed  has  a  real  envelope10,  which  has 
already  been  assumed. 


Adapting  the  discussion  in  Cook  and  Bernfeld,  and  noting  equation 
(2-2),  the  variance  for  the  time-stretch  estimate  is  bounded  as  follows: 


E{(s-s)2} 


“'•o> 

£2h 


(5-7) 


where  h  is  the  input  SNR,  and  £  is  the  RMS  signal  duration,  given  by 
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|  /  t2jf(t)|2dt 
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where  E  is  the  signal  energy,  and  f(t)  is  the  unsampled  signal  envelope. 
Following  the  example  given  in  the  paper  by  Ricker,2  if  the  real 
envelope  has  constant  amplitude,  w  ,  and  bandwidth  B»  the  energy 

O  T 

becomes 
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so  that 
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t  dt 


Et(s-s)2}  > - -  (5-8) 

(o>oT)^h 

If  the  standard  deviation  corresponding  to  (5-8)  is  considered  a 
useful  estimate  of  the  magnitude  of  the  error  in  the  time-stretch 
estimate,  then  (5-4)  becomes 
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At  -  At 


(1  ±  - )  -  At(1  ±  £ - ) 


/h  u  sT 
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so  that  the  bias  term  is  seen  to  be  negligible  for  most  situations,  where 


✓h  «  T  »  /6.  Using  the  standard  deviation  in  (5-6)  yields 


TT 


(5-9) 


This  expression  appears  somewhat  disappointing ,  but  in  practice , 
this  large  bias  estimate  is  often  due  to  the  crudeness  of  the  bound, 
rather  than  to  poor  performance  of  the  estimator.  In  many  situations, 
the  time  difference  of  arrival  is  known  to  have  an  upper  bound  (i.e., 
At<At  ) ,  and  for  a  given  signal ,  one  can  use  this  "worst  case"  time 

SI3X 

difference  in  (5-5)  along  with  the  estimate  of  |s-s|  in  (5-8)  to  calcu¬ 
late  an  estimate  of  the  bias,  .  Indeed,  this  estimator  is  inherently 
limited  to  estimating  delays  that  fall  within  a  range  corresponding  to 
the  time  required  for  the  carrier  to  complete  one  cycle,  i.e. 
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In  essence,  a  priori  knowledge  of  Che  range  of  possible  cime  dif¬ 
ferences  is  required  so  that  one  can  assign  the  proper  time-delay  value 
to  an  observed  phase  difference.  If,  for  example.  At  is  known  to  be 

small  and  positive  then  At  .  will  be  selected  as  zero.  At  «  — 

min  max  f 

o 

and  the  phase,  $,  is  defined  on  (0,2ir). 

If,  for  At  -  At  the  coefficients  f,.  and  f».  are  nearly  the  same 
max  li  Ai 

then  the  imaginary  part  of  the  product  of  the  two  sums  on  the  left-hand 

side  of  (5-5)  will  be  nearly  zero.  For  example,  if  the  carrier  frequency 

is  30kHz,  the  modulation  is  a  850Hz  pure  tone,  the  signal  duration  is 

600  msec,  the  sampling  frequency  is  1700Hz,  and  t  .  *  0,  At  *  , 

min  max  r 

o 

then  the  bias  in  phase  is  calculated  numerically  for  real  envelope 
signals  to  yield  the  following  table: 


TABLE  5-1 

Phase  Bias  for  Real  Envelope  Signal 


Input  SNR 

A%ax  ( vsec ) 

|4>bi  (degrees) 

20  dB 

0.0333 

“  0 

0.333 

9.49  x  10-6 

3.33 

1.49  x  10-4 

33.3  («  l/f0) 

1.44  x  10"3 

10  dB 

0.0333 

3.43  x  10-;> 

0.333 

1.55  x  10-4 

3.33 

H- 
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u> 

33.3 

1.47  x  10-2 

0  dB 

0.0333 

1.97  x  10-14 

0.333 

1.54  x  10“2 

3.33 

1.63  x  JO-2 

33.3 

_ 

1.65  x  10“1 

N\V.v 
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>V'; 


One  would  expect  that  for  more  complicated  real-envelope  waveforms 
with  modulation  frequencies  no  greater  than  850Hz  the  bias  values  would 
be  somewhat  less  than  those  of  Table  5-1 ,  since  for  lower  modulation 
frequencies  the  coefficients  f^  and  f2^  would  be  more  nearly  the  same. 
Hence,  the  simplified  numerical  computations  yielding  Table  5-1  would 
provide  a  means  of  estimating  the  bias  error  of  receivers  utilizing 
more  complicated  real-envelope  waveforms  of  a  given  bandwidth. 

Because  the  method  used  to  generate  Table  5-1  appears  to  be  of  a 
more  general  usefulness,  the  process  is  parameterized,  and  the  FORTRAN 
source  code  is  included  in  Appendix  C.  The  parameterization  is  useful 
because  with  it  one  can  perform  the  computations  once  for  a  whole  class 
of  signals,  rather  than  having  to  repeat  the  computations  for  each 
specific  signal  used.  This  parameterization  is  explained  in  the  comments 
of  the  FORTRAN  source  code. 

If  the  signal  used  is  not  a  real-envelope  signal  as  was  assumed 
above  then  the  above  procedure  is  not  valid,  but  the  bias  can  still  be 
evaluated  numerically  in  the  same  fashion  for  each  specific  signal  using 
the  maximum  time  difference  and  the  standard  deviation  obtained  from  the 
Cramer-Rao  bound.  It  is  mentioned  here  that  if  the  signal  does  not  have 
a  real  envelope,  then  in  general  the  estimates  of  T  and  s  will  in  general 
be  coupled.  In  this  case  the  bound  given  by  (5-7)  will  still  be  valid, 
but  will  not  be  as  tight  as  possible.  A  tighter  bound  that  accounts  for 
x-s  coupling  is  discussed  in  Van  Trees  and  in  Cook  and  Bernfeld. 

It  is  desired  to  combine  the  results  of  this  chapter  with  those  of 
previous  chapters  in  order  to  come  up  with  a  rule  of  thumb  measure  of  the 
total  performance  of  the  estimator.  Here  the  noise  will  be  assumed 


uncorrelated  (as  In  Chapter  3).  The  measure  of  performance  to  be 

2 

derived  will  be  E{($-$q)  },  which  represents  the  expectation  of  the 
squared  deviations  of  the  estimator  from  the  true  value,  <pQ .  The 
previous  discussion  in  this  chapter  implies  that  the  effect  of  the  T-s 
mismatch  can  be  taken  into  account  by  replacing  with  where  <t> 

is  a  zero-mean  random  process  with  variance  determined  from  the 
Cramer-Rao  bound  as  discussed  above.  Without  loss  of  generality,  4>o 
will  be  assumed  zero,  so  that  the  variance  of  Figure  3-3  will  now 
represent  E{($-^)  }.  The  densities  of  Figures  3-1  and  3-2  will  here  be 
represented  as  f(4|$^)  in  accordance  with  the  discussion  following 
Equation  3-10.  Denoting  the  marginal  densities  for  the  bias,  <t>^  and  for 
the  total  estimate,  <p ,  as  f($Q)  and  f < ♦) ,  respectively,  the  density  f(<p) 

9 

can  be  expressed  as 


ir 
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“TT 
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where  the  order  of  integration  has  been  interchanged.  Noting  that 
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♦  -  ^  +  2$$b,  the  inner  integral  becomes 
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In  this  equation,  it  is  noted  that  the  first  integral  is  precisely  the 
2 


variance  a  ,  plotted  in  Figure  3-3 ,  the  second  integral  is  equal  to 
unity,  and  the  third  integral  is  precisely  so  that 


£{(♦“♦  )2}  -  /  [a2  +  d>?]f (4>K)d4> 


-7T 


b'“Tb 


-  o2  +  ob2  .  (5-10) 

Hence,  under  the  assumptions  stated,  one  may  merely  add  the 

2  2 
variance  a  obtained  from  the  Cramer-Rao  bound  to  the  variance  or  from 

D 

Figure  3-3  to  obtain  a  rule-of-thumb  value  for  the  deviations  of  the 
estimator  about  the  true  value.  Note  that  while  this  rule-of-thumb 


was  derived  for  the  uncorrelated  noise  case,  it  is  also  valid  for 
the  correlated  noise  case  if  the  signal-to-noise  ratio  is  high  enough 
that  the  mean  of  the  estimator  is  nearly  equal  to  the  true  value,  <(>  . 


This  concludes  the  discussion  of  the  effects  on  the  estimator  of  a 
t-s  mismatch.  To  summarize,  if  the  time  difference  of  arrival  is  small, 
then  the  only  effect  of  a  t-s  mismatch  is  a  slight  bias  as  in  (5-4)  and 

a  decrease  in  the  effective  signal-to-noise  ratio.  If  the  time 

difference  of  arrival  is  not  small,  then  the  phase  becomes  biased.  A 

crude  estimate  of  this  bias  is  given  in  (5-9).  This  estimate  is 
independent  of  the  properties  of  the  signal.  A  better  estimate  can  be 
obtained  by  assuming  a  maximum  time  difference,  and  by  considering  the 
properties  of  the  signal  employed.  For  real  envelope  signals,  this 
process  can  be  parameterized  and  computed  via  the  program  in  Appendix  C. 
For  more  general  signals,  one  must  use  the  specific  signal  to  calculate 
the  phase  of  (5-5).  A  mismatch  in  t  still  affects  only  the  effective 

sigal-to-noise  ratio. 

2 

The  variance,  o^,  can  be  added  to  the  variances  given  in  the  graphs 
in  the  previous  chapters  to  provide  a  rule-of-thumb  measure  of  the 


estimator  performance 


CHAPTER  6 


SUMMARY  AND  CONCLUSIONS 

It  has  been  shown  in  the  previous  chapters  that  wideband  signals  may 
be  used  for  phase-comparison  time  delay  estimation  provided  the  signals 
employed  have  a  symmetric  power  spectrum.  Chapter  3  analyzed  the 
performance  of  the  estimator  operating  in  uncorrelated  noise,  while 
Chapter  4  discussed  the  performance  in  correlated  noise.  The 
uncorrelated  assumption  is  often  used  in  practice,  where  correlation 
properties  of  the  noise  may  not  be  available. 

It  was  shown  that  in  correlated  noise  the  estimator  becomes  biased, 
but  approaches  an  unbiased  estimator  for  favorable  signal-to-noise 
ratios.  For  signal-to-noise  ratios  in  which  the  densities  derived  are 
valid,  the  variances  fall  off  as  and  the  sample  means  and  variances 

approach  the  derived  theoretical  values.  By  expressing  the  number  of 
input  samples,  N,  as  the  time-bandwidth  product  of  the  signal,  the 
estimator  was  shown  to  use  the  processing  gain  of  large  time-bandwidth 
product  signals  to  reduce  the  estimator  variance. 

In  Chapter  5,  the  effects  of  a  r-s  mismatch  were  examined.  Central 
to  this  is  the  Cramer-Rao  lower  bound  on  the  time-stretch  estimate.  A 
rule-of-thumb  was  given  for  estimating  the  squared  deviations  of  the 
estimator  from  the  true  value. 

This  method  of  time  delay  estimation  has  advantages  over  existing 
methods  in  that  it  can  handle  Doppler  shifted  channels  with  ease,  and 
that  it  can  identify  and  estimate  the  time  delay  for  resolvable  scat- 
terers  each  with  distinct  t-s  values. 
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Further  work  remains  in  evaluating  the  estimator  for  multiple  point 
channels,  and  investigating  possible  Interference  phenomena  between  the 
different  point  reflectors.  Also,  work  remains  in  designing  signals  with 
desirable  characteristics  under  the  new  constraint  that  the  signal  s 
spectrum  be  symmetric. 
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Appendix  A  -  Correlated  Noise  Model 


This  discussion  of  the  noise  correlation  coefficients  is  taken 
directly  from  Merchant**.  The  noise  model  assumed  here  consists  of  two 
components:  an  incoherent  component  plus  a  coherent  component,  i.e.. 


"ki  "  \l  +  nki 


k  -  1,2 


(B-l) 


where  the  superscript  ”i”  denotes  the  incoherent  component,  and  the 
superscript  "c"  denotes  the  coherent  component  of  the  noise.  The  vari¬ 
ances  of  the  real  and  imaginary  parts  of  the  incoherent  and  coherent 
2  2 

components  are  and  a ,  respectively,  where 

2  2  2 
a  *  a.  +  o  . 
n  i  c 


The  incoherent  noise  component  has  the  following  covariance  properties 


E{xli  X2j} 
E{xU  y2j1} 


E{yu  y2ji} 

E{x2i  yl5} 


Pi°i  6iJ 


where  is  the  in-phase  correlation  coefficient  of  the  incoherent 


noise. 


The  coherent  noise  is  assumed  to  te  the  same  in  each  channel  ex- 


br  ■~~s 


para: 


cept  for  a  phase  factor: 


If  the  ratio  of  the  coherent  noise  power  to  the  total  noise  power  is 


so  that  the  noise  correlation  coefficients  used  in  Chapter  4  become 

Etx-x-} 

p  *  — - —  -  r  cos*  +  (l-r)p. 

_  l  c  i 


E{x,y2} 

X  - = —  -  r  sin*  . 

o  2  c 

n 


If  the  incoherent  noise  is  uncorrelated  between  channels  (i.e. ,  P ,  ■ 

P  and  X  can  be  considered  as  the  real  and  Imaginary  components  of  a 

-1* 

"complex  correlation  coefficient,"  y  ■  re  c. 


APPENDIX  B 


Change  in  Correlation  Properties 
Due  to  Complex  Multiplication 

Consider  two  circularly  symmetric  complex  process 

ni  "  *i  +  3*1 
n2  ”  x2  +  iy2 

with  correlations 

Etx^}  -  E{y1y2>  -  po2 

2 

E{x1y2>  •  — E{x2y ^ >  -  Xo  . 

If  these  processes  are  each  multiplied  by  complex  numbers  of  unit 
norm,  it  is  desired  to  find  the  new  correlation  properties.  Let 

Uj  +  jvJ  *  (a  +  jb)(x1  +  jyj)  -  aXj  -  by1  +  .KbXj  +  ay^ 

u2  +  jv2  -  (c  +  jd)(x2  +  jy2)  -  cx2  -  dy2  +  j(dx2  +  cy2) 

where  (a  +  jb)  and  (c  +  jd)  represent  arbitrary  complex  numbers  of  unit 

norm.  Then 

Etu^}  -  E{ (aXj  -  byj)(cx2  -  dy2>} 

-  E{acXjX2  +  bdy^j  -  adx^  -  bcy^} 

■  l(ac  +  bd)p  +  (be  -  ad)X)o2  . 


•.'.V.V.V 
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Similarly, 


E{v^v2)  “  [(ac  +  M)o  +  0>c  -  ad)X}a2 

E{u^v2>  ■  l(ad  -  bc)p  +  (ac  +  bd)X]a2 

^Vlu2^  "  “  ^a<*  ”  ^c)P  +  (ac  +  bd)X]a2 

E(ulvl ^  "  E{u2v2>  -  0  . 

Two  special  cases  are  worthy  of  consideration.  First,  if  both 
processes  n^  and  n2  are  multiplied  by  the  same  number,  then  a  ■  c  and 
b  »  d,  so  that 

E{ujU2>  *  E{v^v2)  *  (a2  +  b2)pc2  ■  pa2 

E (“i v2 )  "  ~E{VjU2}  »  (a2  +  b2)Xo2  -  Xo2 
2  2  2  2 

where  a  +  b  ■  c  +  d  »  1  via  the  assumption  of  unit  norm  multipliers. 
Hence  if  each  process  is  multiplied  by  the  same  complex  number,  the 
correlation  pronertles  do  not  change. 

Another  special  case  is  that  in  which  only  one  channel  is  multiplied 
by  a  complex  number,  e.g.,  a+jb«l(a*l,b*  0).  In  this  case 

E{ujU2>  ■  E{VjV2>  *  (cp  -  dX)o2 
E<ui v2>  ■  EtVjUj}  -  (dp  +  cX)a2  . 

It  is  noted  that  the  general  case  can  always  be  broken  down  into 
consecutive  application  of  the  two  special  cases  (i.e.,  let  c  +  jd  * 


o  n 


APPENDIX  C 


Fortran  Source  Code  for  Phase  Bias  of  Real-Envelope  Signals 

*  This  program  calculates  the  phase  bias  for  real  envelope 

*  signals  as  per  pp.  54-55  of  this  thesis  (Hatlestad,  J.D., 

*  "Phase  Comparison  Time  Delay  Estimation  Using  Wideband 

*  Signals"  The  Pennsylvania  State  University,  December, 

*  1985). 

* 

*  The  input  parameters  are  as  follows: 

* 

*  "time  bandwidth  product"  (FBTEND)  ~  This  parameter  is 

*  the  time  bandwidth  product  of  the  signal.  Since 

*  the  signal  is  assumed  to  have  a  real  envelope  the 

*  bandwidth  used  to  calculate  this  must  be  equal  to 

*  TWICE  the  maximum  modulation  frequency. 

* 

*  "ratio  of  bandwidth  to  carrier  frequency"  (BWRAT)  —  See 

*  at'ove  note  for  "time  bandwidth  product." 

* 

*  "ratio  of  time  difference  to  modulation  period" 

*  (FBDELT)  —  This  parameter  is  the  ratio  of  the 

*  maximum  time  difference  to  be  estimated  to  the 

*  modulation  period  of  the  signal  (1/fmax). 

*  See  page  55  of  thesis. 

* 

*  "ratio  of  sampling  frequency  to  bandwidth"  (FSRAT)  — 

*  Must  be  greater  than  or  equal  to  one. 

* 

*  "input  signal  to  noise  ratio"  (SNR)  —  Signal  to  noise 

*  ratio  in  dB. 

* 

************************************************************ 

COMPLEX  SUMI , SUMJ , DSUM 
PI  -  4 . 0*ATAN (1.0) 

Query  user  for  signal  parameters 
WRITE (6, 10) 

10  FORMAT ( $ , "  Enter  time-bandwidth  product  of  signal:  ") 
READ (5, 15)  FBTEND 
15  FORMAT (F20. 10) 

FBTEND  -  FBTEND/ 2.0 
WRITE(6,20) 

20  F0RMAT($,"  Enter  ratio  of  bandwidth  to  carrier 

&  frequency:  ' ) 

READ (5, 15)  BWRAT 
BWRAT  -  BWRAT/ 2 . 0 
WRITE ( 6,30) 

30  F0RMAT($,~  Enter  ratio  of  time  difference  to  modulation 
&  period  ' ) 

READ (5,15)  FBDELT 
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WRITE (6,40) 

40  FORMAT (  $ , "  Enter  ratio  of  sampling  frequency  to 

&  bandwidth:  ' ) 

READ (5, 15)  PSRAT 
WRITE ( 6,50) 

50  F0RMAT($,'  Enter  input  signal-to-noise  ratio  (dB):  ') 

READ (5,15)  SNR 
C 

C  Convert  SNR  from  dB  to  decimal 
SNR  -  10.0**(  SNR/10.0  ) 

C 

C  Compute  auxiliary  variables  (in  radians) 

NP  -  NINT(  FBTEND*2.0*FSRAT  ) 

FOTEND  -  FBTEND/BWRAT  !  Product  of  carrier  and 

C  signal  duration 

FBTINCR  -  FBTEND/FLOAT(NP)  t  Product  of  bandwidth  and 
C  time  Increment 

FOTINCR  •  FOTEND /FL OAT (NP)  !  Product  of  carrier  and 
C  ti-e  increment 

C 

DELS  -  SQRT ( 6 . 0/SNR) / ( 2 . 0*PI*F0TEND )  !  Expected  error 
C  for  s  estimate 

S  -  1.0  +  DELS 
C 

C  Convert  arguments  of  trigonometric  functions  to  radians 
FBTINCR  -  2 . 0*PI*FBTINCR 
FOTINCR  -  2 . 0*PI*F0TINCR 
FBDELT  -  2.0*PI*FBDELT 
C 

C  Intlalize  sums 

SUMI  -  CMPLX(0. 0,0.0) 

SUMJ  -  CMPLX(0. 0,0.0) 

C 

C  Compute  first  sum 
DO  I-0.NP 

C0SC0SI-C0S(  I*FBTINCR  )*C0S(  S*I*FBTINCR  ) 

-SDMI-SDMI  +  C0SC0SI*CEXP(  CMPLX ( 0 . 0 , DELS*I*F0TINCR)  ) 
ENDD0 
C 

C  Compute  second  sum 
DO  J-0,NP 

C0SC0SJ-C0S(  J*FBTINCR  )*C0S(  S* ( J*FBTINCR+FBDELT)  ) 
SUMJ* SUMJ+C0 SCO SJ*CEXP (  CMPLX ( 0 . 0 , -DELS* J*F0TINCR)  ) 
ENDDO 
C 

C  Multiply  sums 

DSUM  -  SUMI*SUMJ 
C 

C  Compute  phase  bias 

PHIB  -  ATAN(  AIMAG( DSUM) /REAL (DSUM)  ) 

PHIB  -  PHIB*180.0/PI 
WRITE( 6 , 60 )  PHIB 

60  F0RMAT( '  PHIB  (DEGREES)  -  ',E10.4) 

C 

STOP 

END 


irj  rrr 
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